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Abstract 

We  present  a  luetliod  for  solving  the  following  prohlein:  Cliven  a  set  of  data  points  scattered  in 
three  dimensions  and  an  initial  triangnlar  mesh  J/o.  produce  a  mesh  M .  of  the  same  toi)ological 
type  as  A/q.  that  fits  the  data  well  and  has  a  small  number  of  vertices.  Onr  ap])ioach  is  to  minimize 
an  energy  function  that  explicitly  models  the  comireting  dc'sires  of  conciseness  of  lepresentation 
and  fidelity  to  the  data.  Wo  show  that  mesli  o|)timizalion  can  be  effectively  used  in  at  least 
two  applications:  surface  reconstruction  from  unorganizt'd  |)oints.  and  mesh  sim))lification  (the 
reduction  of  the  nutiiber  of  vertices  in  an  initially  den.se  mesh  of  triangles). 

CR  Categories  and  Subject  Descriptors:  l.J.')  [Computer  Craphicsl;  (  ompntational 
Geometry  ami  Object  .Modeling. 

Additional  Keywords:  Geometric  .Modeling.  Surface  I'ittinti.  I'hiee- Dimension. d  Shape  Re¬ 
covery.  Range  Data  .\nalvsis.  .Mothd  Simplification. 


1  Introduction 


The  inesh  oi)tiivi:(ili(ni  itrftblem  coiisiiferc'd  iii  this  paper  can  be  roughly  stated  a^  follo\ts:  Given  a 
collection  of  data  |)oints  .V  in  R*  and  an  initial  triangular  mesh  Mq  near  the  data,  find  a  mesh  M 
of  the  same  topological  type  as  .Un  that  fits  t  h<>  data  woll  and  ha^  a  small  number  (d  vertices. 
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Figure  1;  Examples  of  mesh  optimization.  The  meshes  in  the  top  row  are  the  initial  meshes  Mq; 
the  meshes  in  the  bottom  row  are  the  corresponding  optimized  meshes.  The  first  3  columns  are 
reconstructions;  the  last  2  columns  are  simplifications. 

As  an  example,  Figure  lib  shows  a  set  of  4102  data  points  sampled  from  the  object  shown  in 
Figure  11a.  The  input  to  the  mesh  optimization  algorithm  consists  of  the  points  together  with  the 
initial  mesh  shown  in  Figure  11c.  The  optimized  mesh  is  shown  in  Figure  12c.  Notice  that  the 
sharp  edges  and  corners  indicated  by  the  data  have  been  faithfully  recovered  and  that  the  number 
of  vertices  has  been  significantly  reduced  (from  1572  to  163). 

To  solve  the  mesh  optimization  prolilem  we  minimize  an  energy  function  that  captures  the  com¬ 
peting  desires  of  tight  geometric  fit  and  compact  representation.  The  tradeoff  between  geometric 
fit  and  compact  representation  is  controlled  via  a  u.ser-selectabie  parameter  r,.pp.  A  large  value  of 
Crfp  indicates  that  a  sparse  repre.sentation  is  to  be  strongly  preferred  over  a  dense  one.  usually  at 
the  expense  of  degrading  the  fit. 

We  use  the  input  mesh  Mq  as  a  starting  point  for  a  non-linear  optimization  process.  During 
the  optimization  we  vary  the  number  of  vertices,  their  positions,  and  their  connectivity,  .\lthough 
we  can  give  no  guarantee  of  finding  a  global  minimum,  we  have  run  the  method  on  a  wide  variety 
of  data  sets;  the  method  has  produced  good  results  in  all  cases  (see  Figure  1 ). 

We  see  at  least  two  applications  of  mesh  optimization;  surface  reconstruction  and  mesh  simpli¬ 
fication. 

The  problem  of  surface  reconstruction  from  sampled  data  occurs  in  many  scientific  and  engi¬ 
neering  applications.  In  [2],  we  outlined  a  two  phase  procedure  for  reconstructing  a  surface  from  a 
set  of  unorganized  data  points.  The  goal  of  phase  one  is  to  determine  the  topological  type  of  the 
unknown  surface  and  to  obtain  a  crude  estimate  of  its  geometry.  .\n  algorithm  for  phase  one  was 
described  in  [.5].  The  goal  of  phase  two  is  to  improve  the  fit  and  reduce  the  number  of  faces.  Mesh 
optimization  can  be  used  for  this  purpose. 

Although  we  were  originally  led  to  consider  the  mesh  optimization  problem  by  our  research 


on  surface  reconstnictiou.  the  algorithm  we  liave  (levelo|)<'cl  can  also  he  applied  to  the  problem  of 
mesh  simplification.  .Mesh  simplilication.  as  considered  by  Tiirk  [lb]  and  Schroeder  et  al.  [11],  refers 
to  the  problem  of  reducing  the  number  of  faces  in  a  dense  mesh  while  minimally  perturbing  the 
shape.  Mesh  optimization  can  be  vised  to  solve  t  his  problem  as  follows;  sample  data  points  X  from 
the  initial  mesh  and  tise  the  initial  mesh  as  the  starting  point  Mo  of  the  optimization  procedure. 
For  instance.  Figure  14a  shows  a  triangular  approximation  of  a  minimal  surface  with  2032  vertices. 
.Application  of  our  mesh  optimization  algorithm  to  a  sample  of  6752  points  (Figure  14b)  from  this 
mesh  produces  the  meshes  shown  in  Figures  14c  (487  vertices)  and  14d  (239  vertices).  The  mesh 
of  Figure  14c  corresponds  to  a  relatively  small  value  of  Crep-,  and  therefore  has  more  vertices  than 
the  mesh  of  Figure  14d  which  corresponds  to  a  somewhat  larger  value  of  Cr^p- 

The  principal  contributions  of  this  paper  are: 

•  It  presents  an  algorithm  for  fitting  a  mesh  of  arbitrary  topological  type  to  a  set  of  data  points 
(as  opposed  to  volume  data.  etc.).  During  the  fitting  |)rocess.  the  numiver  ttiid  connectivity 
of  the  vertices,  as  well  as  their  positions,  are  allvwed  to  vary. 

•  It  casts  mesli  simplification  as  an  o|)timization  problem  with  an  energy  function  that  directly 
measures  deviation  of  the  litial  mesh  from  the  original.  .\s  a  consetiuence.  the  final  mesh 
naturally  adapts  to  curvature  variations  in  the  original  mesh. 

•  It  demonstrates  how  the  algorithm's  ability  to  recover  sharp  edges  and  corners  can  be  ex¬ 
ploited  to  atitomatically  segtnent  the  final  mesh  into  smooth  connected  comitonents  (see 
Figure  Fie). 

2  Mesh  Representation 

Intuitively,  a  ;;)c-.s7(  is  a  piecewise  linear  surface,  consisting  of  triangu'ar  faces  pastt'd  together  along 
their  edges.  For  our  ivnrposes  it  is  important  to  maintain  the  distinction  between  the  connectivity 
of  the  vertices  and  their  geometric  positions.  Formally,  a  mesh  M  is  a  pair  (  A'.D.  where:  A'  is  a 
siinplicial  coinpifx  representing  the  ( onnectivity  of  the  vertices,  v'dges.  and  faces,  thus  determining 

the  topological  type  of  the  mesh:  1’  =  [v] . v,„}.  v,  €  R*  is  a  set  of  vertex  positions  defining 

the  shape  of  the  mesh  in  R'  (its  geometric  realization). 

.A  simplicial  com))lex  K  consists  of  a  s(>t  of  vertices  {1 . in},  together  wit  h  a  set  of  non-empty 

subsets  of  the  vertices,  called  the  simplices  of  A',  such  that  any  set  consisting  of  exactly  oite  vertex 
is  a  simplex  in  A'.  ;ind  ew'rv  non-v'inpiy  sub.set  of  a  simplex  in  A'  is  again  a  sinijrlex  in  A’  (cf. 
Spanier  [15]).  rite  O-simplices  {/}  €  A'  ;ire  called  vertices,  the  l-sim|)lices  {i.j}  t  A’  are  called 
edges,  and  the  2-sim()lices  {l.j.h}  €  A’  tire  called  faces. 

•A  geometric  realization  of  a  mesh  as  a  surface  in  R*  can  be  obtained  as  follows.  For  :v  given  sim¬ 
plicial  complex  A’,  loim  its  lopoloiiiml  r<nlizalii>n  |A'|  in  R'"  l)y  identifying  the  vertices  {1 . ;??} 

with  the  standard  l)asis  vc'ctors  {ej . e  }  of  R'".  For  <'ach  -imi>lex  >  •£  A'  h'l  f-i  denote  i  li(> 

convex  hull  of  its  vertices  in  R'".  and  let  |A’|  =  1-^1.  l.t't  o  :  R'"  —  R-*  lx-  the  linear  map  that 

sends  the  /-th  standard  basis  \ector  e,  R"'  to  V,  fc  R  (see  Figure  2). 


Simplicial  complex  A' 
vertices;!  1 },  {2},  {3} 
edges;  {1,2},  {2,  3} ,  { 1,  3} 
faces;  { 1 ,  2,  3 } 

Topological  realization  \K\  Geometric  realization 


Figure  2:  Example  of  mesh  represetitatioii:  a.  mesh  consisting  of  a  single  face. 


The  geometric  veolizalioo  of  M  is  f  lie  image  oi-(  |  A’l).  where  we  write  the  map  as  oy  to  emphasize 

that  it  is  fully  specified  by  the  set  of  vertex  positions  \’  -  {vi . v„,}.  The  map  oy  is  called 

an  embedding  if  it  is  l-l.  that  is  if  c;>r(|A’j)  is  not  self-intersecting.  Only  a  restricted  set  of  vertex 
positions  I-'  result  in  oy  being  an  enibedding. 

If  0v  is  an  embedding,  any  point  p  €  oi-(|A'l)  can  be  parameterized  by  finding  its  unique 
pre-image  on  jA’].  The  vector  b  €  |  A'|  with  p  =  oi'(b)  is  called  the  bnrycentric  coordinate  rector  of 
p  (with  re.spect  to  the  simplicial  (omplex  A').  .N'ole  that  barycentric  coordinate  vectors  are  convex 
combinations  of  standard  basis  vectors  e,  6  R'"  corresponding  to  the  vertices  of  a  face  of  A  .  .Any 
barycentric  coordinate  vector  has  at  most  three  non-zero  I'ntries:  it  has  onl>'  two  non-zero  entries 
if  it  lies  on  an  edge  of  1A'|.  and  only  one  if  it  is  a  vertex. 


3  Definition  of  the  Energy  Function 

Recall  that  the  goal  of  mesh  o|)timizat ion  is  to  obtain  a  mesh  that  inovides  a  good  fit  to  the  point 
set  X  and  has  a  small  number  of  vertices.  We  find  a  simplicial  complex  h  and  a  set  of  vertex 
positions  \  defining  a  mesh  M  =  (  A  .  \  )  tliat  minimizes  t  Ih>  <'nergy  function 

A(  A  .  1  )  =  A./i.s((  A  .  1  )  -f-  A,.fp(  A  )  un/l  A  .  1  ). 

The  first  two  terms  correspond  to  the  two  stated  goals:  the  third  term  is  motivated  below. 

The  distance  energy  is  eipial  to  the  sum  of  sipiared  distances  from  the  points  .\  = 

(xi . x„}  to  the  mesh. 

A.n.w(  A'.  T)  =  ^d'‘(x,.oi-(|A'l)). 


t 


The  represeiUa.lioii  cMiersy  Ertp  penalizes  ineslies  with  a  lai»;o  iiuinlter  ol  vertices.  It  is  set  to 
be  proportional  to  the  miinher  of  vertices  rn  of  K: 

The  optimization  allows  vertices  to  l)e  both  added  to  and  removed  from  the  mesh.  W  hen  a  vertex 
is  added,  the  distance  energy  is  likely  to  be  reduced:  the  term  Erep  makes  this  operation  incur 
a  penalty  so  that  vertices  are  not  added  indelinitely.  Similarly,  one  wants  to  remove  \  ertices  from 
a  dense  mesh  even  if  increases  slightly:  in  this  case  E^ep  acts  to  encourage  the  vertex  removal. 
The  user-specified  parameter  Crep  provides  a  controllable  trade-off  between  fidelity  of  geometric  fit 
and  parsimony  of  representation. 

We  discovered,  as  others  have  l)efore  us  [7],  that  minimizing  E.i,st  +  Er^p  does  not  produce  the 
desired  results.  .\s  an  illustration  of  what  can  go  wrong.  Figure  1  Id  shows  the  result  cri  minimizing 
E.jist  alone.  The  estimated  surface  has  sevi'ial  spikes  in  regions  where  there  is  no  data.  These  spikes 
are  a  manifestation  of  the  fundamental  problem  that  a  minimum  of  E.nst  +  Ertp  may  not  exist. 

To  guarantee  t  he  existence  of  a  minimum  (see  .Vppendi.x  .\ .  1 ).  we  add  t  he  t  hire!  term,  t  he  s])ring 
energy  £',,prtn-;-  It  places  on  each  edge  of  the  mesh  a  spring  of  rest  length  zero  and  s|)ring  constant 
k: 

Espr,n:,(  l\.\')=  Yi  ''I!'''  " 

It  is  worthwhile  emphasizittg  that  the  spring  energy  is  not  a  smoothness  ])enalty.  Our  intent 
is  not  to  penalize  sharj)  dihedral  angles  in  the  mesh,  since*  such  features  may  be  i)rcsent  in  the 
underlying  surface  and  should  be  recove'red.  We  vi(>w  E^p,-,,,,,  as  a  i('gnlarizing  te'rm  that  helps  guide 
the  optimization  to  a  elesirable  local  minimum.  .\s  the  o()t imizat ion  converges  to  tin*  solution,  the 
magnitude  of  £,,,,,.1-17  <’an  be  gradually  reduced.  We  return  to  this  issue  in  Section  1.1. 

For  some  ap|)licat  ions  we*  want  t  he  |)roce<lnre  to  be  sea  k'-in  variant  .which  is  ('([ui  vak'iit  t  o  (k'fining 
a  nnitless  energy  funct  ion  A  .  To  achieve  invariance  under  Fnclidean  mot  ion  and  uniform  scaling,  t  he 
jjoints  X  and  the  initial  mesh  .f/o  are  |)re-scaled  uniformly  to  (it  in  a  unit  cube,  .\fter  optimization, 
a  post-processing  step  can  ntido  this  initial  transformation. 

4  Minimization  of  the  Energy  Function 

Our  goal  is  to  minimize  the  energy  function 

A(  A  .  1  )  =  /■-./,.•(/{  A  .  I  )  -h  trfp(  A  )  +  A.spriiivl  A  .  \  ) 

over  the  set  X  of  sim|)licial  complex('s  A’  hotneomorphic  to  the  initial  simplicial  complex  A'q.  and 
the  vertex  positions  \’  defitiing  the  embe<lding.  We  now  present  an  outline  of  our  optimization 
algorithm,  a  pseudo-code  version  of  w  hich  appears  in  Fignia*  d.  The  details  are  deferrc'd  to  the  next 
I  wo  subsections. 

To  minimize  A'(  A'.  \ ' )  over  both  A  and  F.  we  part  it  ion  t  he  problem  into  t  wo  nested  subproblems: 
an  inner  minimization  over  F  for  li.xed  simplicial  complex  A  .  and  a  outer  minimization  over  A’. 


In  Section  4.1  wo  doscribe  an  alj^oritlun  liiat  solves  t  lie  inner  minimization  |)rol)leni.  It  finds 
E(I\)  =  mini’  E(K.  1).  the  energy  of  the  best  jrossilrle  <‘mbedding  of  the  fixed  simplicial  complex 
K .  and  the  corresponding  vertex  irositions  I',  given  an  initial  guess  for  V'.  This  corresjx.nds  to  the 
procedure  OptimizeVertexPositions  in  Figtire  4. 

Whereas  the  inner  minimization  is  a  continuous  optimization  prof  em.  the  outer  minimization 
of  E{K)  over  the  simplicial  complexes  K  6  A.’  (procedure  OptimizeMesh)  is  a  discrete  optimization 
problem.  An  algorithm  for  its  .solution  is  presented  in  Section  1.2. 

The  energy  function  E{K.V]  depends  on  two  parameters  c,.ep  and  k.  The  parameter  Crep 
controls  the  tradeoff  between  conciseness  and  fidelity  to  the  data  and  should  be  set  by  the  user. 

The  parameter  k.  oii  the  other  hand,  is  a  regularizing  parameter  that,  ideally,  would  be  chosen 
automatically.  Our  method  of  setting/;  is  described  in  Section  1.4. 

4.1  Optimization  for  Fixed  Simplicial  Complex 
( Procedure  OptimizeVertexPositions ) 

In  this  section,  wo  consider  the  proltlem  of  finding  a  set  of  vertex  positions  1  tliat  minimizes  the 
energy  function  E{h'.\')  for  a  given  simplicial  complex  l\ .  .\s  7:,fp(  A  )  does  not  deitend  on  V'.  this 
amounts  to  minimizing  d  A'.  1')  +  E, E A'). 

To  evaluate  the  distance  energy  A./.a/l  A  .  V).  it  is  nece.s.sary  to  compute  the  distance  of  each 
data  point  x,  to  .1/  =  oi-((A'|).  Fach  of  these  distances  is  itself  the  solution  to  the  minimization 
problem 

r/-(x,.or(|A'l))  =  min  ||x,  -  or(b,)||‘‘. 

in  which  the  unknown  is  the  l/aryceiit tie  coordinate  vector  b,  G  |A’l  C  R’"  of  the  projection  of  x, 
onto  M  (Figure  I).  1‘hns.  minimizing  /-.(A.  1)  for  fixed  A’  is  ('(|uivaleut  to  mittimiziug  the  new 
ob  jective  function 

n 

/■.'(  A'.  WB)  =  11^'  “  'A  l  b,  llj-  +  A'.  I ' ) 

/=! 

/} 

'=i  l.,.AleK' 

over  the  vertex  positions  I '  =  {vj . v,„}.v,  G  R*  and  t  he  l/arycent  ric  coordimites  B  =  (bj . bT,}.b,G 

lA'IcR'T 

To  solve  this  optimization  proltlem  (procedure  OptimizeVertexPositions).  our  tnethod  alternates 
between  two  stibprol/lems: 

1.  For  fixed  vortex  i/ositions  I',  find  optimal  barycentric  coordinate  vectors  B  by  projection 
(procedure  ProjeetPoints). 

2.  For  (i.xed  barycentric  coordinate  vectors  B.  (ind  optimal  \erlex  positions  1'  by  solving  a  linear 
least  squares  irroblem  (procedure  ImproveVertexPositions). 


(i 


OptimizeMesh(  A’o-^  o)  { 

K  ;=  A'o 

V'  :=  OptimizeVertexPositionsi  A'o.^  o) 

-  Solve  the  outer  niininiizritioii  probleiu. 

repeat  { 

(A’M  ')  ;=  GenerateLegalMove(  A  ) 

=  OptimizeVertexPositions(  A  M  ') 

if  A(A''.r')  <  /'.'(A.r)  tlieii 
(AW)  :=  (  A'.\-') 
endif 

}  until  converjiietK'o 
return  (A’.V) 

} 

-  Solve  the  inner  optiniizution  jHohleni 

£(  A')  =  ininr  A'(  A‘.  T) 

-  for  fixed  siniplieiul  ronipler  A’. 

OptimizeVertexPositions(  AM  )  { 

repeat  { 

C  onipiite  htn  jp  e  ntrir  cnorfliinjli  s  hi/  jn'itji  rliiiii. 

B  :=  ProjectPointsl  A  .\  ) 

Miniini:(  l:{  h  .\  .  B)  nn  r  1  us/iif/  ronjiif/iiti  i/riidn  n/.^. 

I  ;=  ImproveVertexPositionsl  A’.//) 

}  until  convor,i^en(  e 
return  1 ' 

) 

GenerateLegalMovel  A'.r )  { 

Select  a  legal  move  A'  =>  K' . 

Locally  modify  \  to  ()l)tain  I  '  appropriate  k)r  A  ', 
return  (A'M'') 

} 

Figure  d;  An  ideali/c'd  ps('udo-code  vrusion  of  t  lie  minimization  algorithm. 


Fimiro  1;  l)islniuoor;i  |)oiiii  x,  (Vom  the  mesh. 


Because  we  find  optimal  .solutions  to  hot  h  of  t  li<‘s<'  snl)i>iol)lems.  E(  A  .  \  .  B )  can  never  increase, 
and  since  it  is  l)ounded  from  l)elow.  it  must  converge.  In  principle,  one  could  iterate  until  some 
formal  convergence  criterion  is  met.  Instea<l.  as  is  common.  w(>  perform  a  fi.ved  number  of  iterations. 
.\s  an  example.  Figure  1  le  shows  the  result  of  optimizing  the  mesli  of  Figure  11c  over  the  vertex 
positions  while  holding  the  simjtlicial  complex  fixed. 

It  is  conceivable  that  |)rocedur(‘  OptimizeVertexPositions  returns  a  set  of \ertices  for  which  the 
mesh  is  .self-intersecting,  i.e.  oi-  is  not  an  embedding.  W  hile  it  is  possibh'  to  check  n  posteri¬ 
ori  whether  Of  is  an  emlx'dding.  constraining  the  optimization  to  always  itroduce  an  embedding 
appears  to  be  difficult.  This  has  not  presentt'd  a  probhmi  in  the  (>xamples  we  hiiv('  run. 

4.1.1  Projection  Subproblein 
(Procedure'  ProjeetPoints ) 

File  problem  of  optimizing  /■,(  A.l  ,  U)  over  the  barycentric  coe)rdin.ite  vectors  B  —  {bj . b,,}. 

while  holding  the  vertex  positions  1'  =  {V] . v,„}  ami  the  simplic'd  complex  A’  constant,  de¬ 

composes  into  n  soi)arate  optimization  problems; 

b,  =  argmin  ||x,  -  (b)|| 

helh'l 

In  other  words,  b,  is  the  l)arycenl  tie  coordinate  ex'cior  corresponding  to  tlie  point  p  e  oi(  |  A’| ) 
clo.sest  to  X,. 

A  naive  approach  to  computing  b,  is  to  project  x,  onto  all  of  the  faces  of  .\/.  and  then  find  the 
projection  with  tninimal  distance,  'lu  speed  up  the  projection,  we  first  enter  the  faces  of  the  mesh 
into  a  spatial  partitioning  data  structure  (similar  to  the  one  used  in  [17]).  Then  for  each  point  x, 
only  a  nearby  subsc't  of  the  faces  needs  to  be  considered,  and  the  projection  step  takes  expected 
time  0{n).  lor  additional  speedu|)  we  exploit  coherence  Ix'tween  iterations.  Insli'ad  of  projecting 
each  point  globally  onto  the  mesh,  we  assume  that  a  iioint's  i)rojeclion  lies  in  a  ix’ighborhood  of 
its  projection  in  the  previous  iteration.  Specifically,  we  project  the  ix)int  onto  all  faces  that  share 


•S 


a  vertex  witli  the  previous  face.  Alt lioiij^li  this  is  a  heuristic  that  can  tail,  it  lias  perlonned  well  in 
Itractice. 

4.1.2  Linear  Least  Squares  Subproblem 
( Procedure  ImproveVertexPositions ) 

Minimizing  E(  K .W  U)  over  the  vertex  positions  I'  while  holding  B  and  A  fixed  is  a  linear  least 
squares  problem.  It  decotnposes  into  three  independent  subproblems,  one  for  each  of  the  three 
coordinates  of  the  vertex  positions.  W’e  will  write  down  the  problem  for  the  lirst  coordinate. 

Let  e  be  the  nutnber  of  edges  (  1 -siniplices)  in  A’:  note  tnat  r  is  Oft)))-  Let  v‘  he  the  ni- 
vector  whose  i-th  element  is  the  first  coordinate  of  v,.  l  et  d'  he  the  ( n  +  f  )-vector  whose  lirst  n 
elements  are  the  first  coordinates  of  the  data  points  x,.  and  whose  last  f  elements  are  zero.  With 
these  definitions  we  can  express  the  U'ast  s(piares  |)rohlem  for  the  first  coordinate  as  minimizing 
j|.lv*  -  d'|l‘  over  v’.  I'he  design  matrix  .1  is  an  (n  +  ^ )  m  matrix  of  scalars.  The  first  n  rows  of 
.1  are  the  barycetitric  coordinate  vi'ctors  b,.  i'.ach  of  tin*  trailing  (  rows  contains  2  non-zero  entries 
with  values  \/h  and  — \/k  in  the  columns  corresponding  to  the  indices  of  the  edge's  endpoints. 
The  first  n  rows  of  the  least  sipiares  problem  corrosiiond  to  A’.  1  ’ ).  whih'  the  last  c  rows 

correspond  to  A’.  T).  .\n  important  feature  of  t  he  matrix  .1  is  that  it  contains  at  most  .1 

non-zero  entries  in  each  row.  for  a  total  of  (){ii  +  in)  non-zero  entries. 

To  solve  the  least  siptares  |)rohlem.  w(>  use  the  conjugate  gradient  method  (cf.  [.'1]).  This  is  an 
iterative  method  guaranteed  to  find  the  exiu  t  solution  in  as  many  itmations  as  then'  are  distinct 
singular  values  of  .1.  i.e.  in  at  itiost  :n  iterations.  Tsually  far  fewer  iterations  are  re(|uired  to  get 
a  resitlt  with  acceptable  precision,  for  I'xamph'.  w('  find  that  for  in  as  large  as  If)',  as  few  as  200 
iterations  are  sufficient. 

The  two  time-consuming  operations  in  <‘acli  iteration  of  the  (on jugate  tiiadient  alirorithm  are 
the  multiiilication  of  .1  by  an  [ii  -|- <  i  vi-ctor  and  tin*  multiplication  of  ,1'  by  an  m-\'ector.  Because 
.1  is  sparse,  tliesi'  two  o|)erations  can  lie  execut<>d  iu  (){ii  -f  ml  time.  We  vior<‘  1  in  a  sparsi'  form 
that  requires  only  ()(ii  -|-  m|  space.  I'lius.  an  acceptable  >olution  to  'he  least  s(|uares  iirobletn 
is  obtained  in  0(ii  in)  time.  In  contrast,  a  typical  noniterative  im'thod  for  solving  dense  least 
sqitares  jirobletns.  such  as  (JH  decomiiositiou.  would  re<|uir('  (A("  +  ni)nr  )  time  to  find  an  exact 
sohttion. 

4.2  Optimization  over  Simplicial  Complexes 

(Procedure  OptimizeMesh ) 

To  solve  the  outer  minimization  problem,  miiiimiziim  /-.(A')  over  A’,  we  (h'line  a  m-i  of  three  ele¬ 
mentary  transformations.  irA/r  co/A/p.sf.  itlii<  siilil.  and  <r///(  Mcup.  taking  a  simidiiial  cotniilex  A’ 
to  another  simplicial  complex  A"'  (see  figun'  b). 

\Ve  define  a  hiinl  nini-i  to  In-  the  pplication  of  on<'  of  thi'si'  ('leiiH'ntary  transformations  to  an 
e<lge  of  A  that  leavi's  tin'  topological  type  of  A  unchan»('d.  Ihe  set  of  elementarx  t  laiislormations 
is  complete  in  the  sense  that  mii/  simplicial  complex  in  K  can  be  obtained  from  A'n  through  a 


sequence  of  legal  moves'. 

Our  goal  then  is  to  find  sncIi  a  se(|uenre  taking  us  from  A'o  to  a  minimum  of  La  I\  i.  We  do  this 
using  a  variant  of  random  desccmt :  we  randomly  select  a  legal  move.  K  A’'.  If  A(  A''  >  <  E{  A  ),  we 
accept  the  move,  otherwise  w<‘  try  again.  If  a  large  nnmlx'r  of  trials  fails  to  prodme  an  acceptable 
move,  we  terminate  the  search. 

.More  elaborate  selection  st  rategies.  such  as  steepest  descent  or  simulated  aunealimr.  are  possible. 
.Vs  we  have  obtaim'd  good  results  witli  the  simple  strategy  of  random  desctuif.  we  have  not  yet 
implemented  the  other  strategies. 

Identifying  Legal  Moves  .\n  edge  "plii  t  rau>fortnat  ion  always  a  legal  mo\e.  a>  it  t  an  never 
change  the  topohjgical  ty|)e  of  A',  fin'  oi  lior  two  t  ratisforniatioiis.  on  tlie  other  hand.  c;ui  cause  a 
change  of  topological  type,  so  tests  must  he  iterfortned  to  deleiinine  if  they  are  legtd  mo\es. 

We  define  an  edge  (A  /}  t  A'  to  he  a  hotiiidai ./  <(!(/<  ifit  In  a  subset  of  only  one  fate  {i.j.L}  €  A’, 
and  a  verte.x  {/}  to  he  a  hoimddrn  irrt(.r  if  there  e.xists  a  hot  ndtiry  t'dge  {i.jj  €  A  . 

-■Vn  edge  collap.se  transformation  A'  =>  A'  that  collapses  the  edge  {i.J}  €  A’  is  a  legal  move  if 
and  only  if  the  following  conditions  are  satisfied  (proof  in  .\ppendi.\  .\.2): 

•  For  all  vt'i'tict's  {/.}  adjiiceni  to  ht>ih  •jM  and  \j}  ((/./.}  A  itnd  fj./'}  A  I.  {i-j.k}  is  a 
face  of  A'. 

•  If  {/}  and  \j}  are  both  hoiiiid.try  \'Mlic('s.  {/.,/}  is  ;t  boundary  etige. 

•  A’  has  more  ilian  I  \'eriices  if  iiciiher  (/}  nor  {7}  are  houndar\-  vtutices.  or  A'  h;i>  more  than 
d  vertices  if  either  {/}  or  are  boundary  vertict's. 

'In  fact,  we  prove  in  Appendix  .\.!  that  edge  collapse  and  edge  s|)lii  are  snilicieni:  we  incliule  edge  swap  to  allow 
the  optimization  procedure  to  limner'  lliroiigli  -..1131)  liills  jii  the  energy  function. 


10 


An  edge  swap  t ransibniiat ion  l\  l\'  thal  replacc's  I  lie  edge  [i.j]  €  A  with  {/<'./}  €  A'  is  a 
legal  move  if  and  only  if  {A’./}  ^  A’. 

4.3  Exploiting  Locality 

The  idealized  algorithm  described  so  lar  is  loo  inelficieni  to  he  of  practical  use.  In  this  section,  we 
describe  some  heuristics  which  dramat i<ally  rerliice  the  running  time.  These  heuristics  capitalize 
on  the  fact  that  a  local  change  in  the  strnctiire  of  the  mesli  leaves  the  optimal  positions  of  distant 
vertices  essentially  unchanged. 

4.3.1  Heuristics  for  Evaluating  the  Effect  of  Legal  Moves 

Our  strategy  for  selecting  legal  move's  reepiires  evaluation  of  E(  A  ')  =  min\-  A  (  A"'.  1)  for  a  simplicial 
complex  A"'  obtained  from  A'  through  ;i  legal  move.  Ideally,  we  wotild  use  procedure  OptimizeVer- 
texPositions  of  Section  l.l  for  this  purpo.se.  as  indicated  in  Figure'  3.  In  practice,  however,  this  is 
too  slow.  Insteael.  we  use'  fast  loeal  lieuristics  to  e'stimate  the  effect  of  a  legal  move  on  the  energy 
function. 

Each  of  the  heuristics  is  Itaseel  e)n  ext racting  a  submesh  in  the  neighborhood  of  the  transforma¬ 
tion.  along  with  the  sultse't  eif  the'  elala  pennis  projecting  onto  the'  submesh.  The  change  in  overall 
e'liergy  is  estimated  by  only  consielering  the  contribution  of  the  subteiesh  and  the  corresponding 
point  set.  This  estimate  is  always  pc'ssimist ic.  as  full  optimization  would  only  further  reduce  the 
energy.  Therefore,  the  heuristics  never  suggest  changes  that  will  increa.se  the  true  energy  of  the 
mesh. 


Definition  of  neighborhood.s  in  a  simplicial  complex  To  re-fer  ic)  neisihborhoods  in  a  .simpli¬ 
cial  complex,  we  nec'd  to  iutreiduee'  Mime  further  notation.  We  write'  to  denote  that  simplex 

''  is  a  non-empty  sulise't  or>imple.v  s.  for  simplex  s  6  A'  (I'igure'  li): 

starf.s:  A')  =  G  A'  ;  <  •"'} 

starf.s;  A  )  =  {.s^  G  A  :  3/  G  siarf  '':  A  )  ; 
linkfs:  A' )  =  star!  .s;  A’ )  \  st  ar(  .s;  A’  )• 


Evaluation  of  Edge  Collapse  lo  evaluate  a  transformation  h'  =>  A'  collapsing  an  edge  {i.j} 
into  a  single  ve'rtex  {//}  (Figure'  a),  we'  take  the  subme'sh  to  be'  si  ar(  { / } :  A' )  U  st  ar(  {j } :  A' ).  and 
optimize  over  the  single  verte-x  position  v/,  while  holding  all  othe'r  \e'rte'x  positions  constant. 

Because  we  [X'lforin  only  a  siiiall  iiiiiiiberof  iieralions  (for  re'tisons  of  e'fficie'iicy  ).  the  initial  choice 
of  V/,  greatly  influence's  I  lie  acciirac\'  of  t  lie  result .  riu'rt'lore'.  \\i'  ;it  ti'iitjit  t  liree  o|)timi/.ations.  with 
V/i  starting  at  v,.  v,.  and  -vf  v,  -f-  v,).  and  .iccept  the  be'sl  eme. 
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star{s,K} 


star{s,K} 


/  /  y> 

t  star(t,K)  star{t,K] 

Fisiii'o  <i;  Xojghborliood  siil>sets  ol'  A'. 


7;  Two  loc.il  lo  .su;ii) 


The  edge  collapse  should  ho  allowed  ouly  if  the  now  luosh  does  not  iuuusoct  ilNolf.  Chorkiug  for 
this  would  be  costly:  instead  we  s<'ttte  for  a  less  expensive  li('iiiist ic  check.  If.  after  the  local  opti¬ 
mization.  the  maximum  dihedral  angle  of  the  edges  in  siar({/)}:  A')  is  greater  than  some  threshold, 
the  edge  collapse  is  rejected. 

Evaluation  of  Edge  Split  The  procedure  is  the  same  ;is  for  (ulge  collapse.  <'xcept  that  the 
submesh  is  defined  to  Ite  star!  {/.,/}:  A  ),  and  the  initial  value  of  the  new  vertex  v;,  is  chosen  to  be 

7(v,  +  Vj). 

Evaluation  of  Edge  Swap  i'o  (waluate  an  edge  swa|)  t  r:uisform;>tion  A’  =>  A  '  that  replaces 
an  edge  {i.j}  €  A'  with  {/•'./]  €  A',  we  consider  two  local  optimizations,  one  with  submesh 
star({/.’}:  A'),  varying  vautex  vj..  and  one  with  snbmesh  star({/}:  A'),  varyitig  vertex  v/  (Figure  7). 
The  change  iti  energy  is  taken  to  best  of  these.  .\s  is  the  case  in  evaluatitig  an  erlge  collapse,  we 
reject  the  transformation  if  the  maxitmim  dihedral  angle  after  tlie  local  optimization  exceeds  a 


threshold. 


4.3.2  Legal  Move  Selection  Strategy 
(Procedure  GenerateLegalMove) 

The  simple  strategy  for  selecting  legal  moves  described  in  Section  1.2  can  be  improved  by  e.xploiting 
locality.  Instead  of  selecting  edges  completely  at  random,  edges  are  selected  from  a  candidate  set. 
This  candidate  set  consists  of  all  edges  that  may  lead  to  benelicial  moves,  and  initially  contains  all 
edges. 

To  generate  a  legal  move,  we  randondy  remove  an  edge  from  the  candidate  sc't.  We  first 
consider  collapsing  the  edge,  accepting  the  move  if  it  is  legal  and  reduces  the  total  energy.  If  the 
edge  collapse  is  not  accepted,  ue  then  consider  (>dge  swap  and  c'dge  split  in  that  order.  If  one  of 
the  transformations  is  accc'pted.  we  update  the  candidatt*  set  by  adding  all  neighboring  edges.  The 
candidate  set  becomes  vtuy  iiselul  toward  the  (uid  of  opiimixation.  when  the  fraction  of  beneficial 
moves  diminishes. 

4.4  Setting  of  the  Spring  Constant 

We  view  the  spring  energy  as  a  regularizing  term  that  helps  guide  the  optitnization  process 

to  a  good  minimutti.  The  spring  constant  k  determines  the  contribiition  of  this  term  to  the  total 
energy.  We  have  obtained  good  results  l)y  making  successive  calls  to  procedure  OptimizeMesh.  each 
with  a  different  value  of  «.  according  to  a  schedule  that  gradually  decreases  k. 

As  an  e.xample.  to  obtain  tin*  final  mesh  in  l'igur(>  12c  starting  from  tin'  nn'sh  in  Figure  11c. 
we  successively  set  to  10"'^.  10“  *.  10“  '.  and  I0“''‘  (s('e  Figuia's  1  If- 12c).  This  same  schedule  was 
used  in  all  the  e.xamples. 


5  Results 

5.1  Surface  Reconstruction 

from  the  set  of  poittts  shown  in  Figure  1  lb.  phase  one  of  our  reconst  ntct  ion  algori'  hm  [■')]  produces 
the  mesh  shown  in  1  ignre  11c:  this  mesh  has  the  corn'ct  topological  tyi)e.  but  it  is  rather  dense,  is 
far  away  from  the  data,  and  lacks  the  shar])  featnii's  of  the  origitial  model  (Figiire  1  la).  Fsing  this 
tnesh  as  a  starting  point,  mesh  o|)timization  produces  the  mesh  itt  Figure  12c. 

Figures  13a-FJI  show  two  e.xamples  of  surface  n'coiistriution  from  actual  laser  range  data  (cour¬ 
tesy  of  Technical  .\rt.s.  lledmond.  W.\).  Figures  Ida  and  Fib  show  sets  of  points  obtained  by  sam¬ 
pling  two  physical  ob  jects  (a  distributor  ca])  and  a  golf  club  head)  with  a  laser  range  finder.  The 
outputs  of  phase  oik'  are  shown  in  Figures  l.'lc  and  Idd.  Fhe  holes  present  in  the  surface  of  Fig¬ 
ure  13c  are  artifacts  ol  t  he  data,  as  sell-shadowing  prexu'iited  some  regions  of  the  surface  from  being 
scanned.  .\daptiv('  selection  oi  scanning  paths  pr<'v«'nting  such  shadowing  is  an  ititeresting  area  of 


future  researcli.  lii  this  <;iso.  \\('  iiiamially  filled  the  hoh-s.  It'aviiig  a  single  l)oiuular\’  at  tlie  l)Ottoni. 
Figures  13e  and  idfshow  the  optimized  meshes  ohtaiiied  with  our  algorithm. 

5.2  Mesh  Simplification 

For  mesh  simplification,  we  first  sample  a  set  of  points  randomly  from  the  original  mesh  using 
uniform  random  sam])ling  over  area.  N'ext.  we  add  the  vertices  of  the  mesh  to  this  point  set. 
Finally,  to  more  faitlifully  preserve  the  boundaries  of  the  mesh,  we  sample  additional  points  from 
boundary  edges. 

.\s  an  example  of  mesli  simplification,  we  start  with  t  he  mesh  containing  2032  vertices  shown  in 
Figure  I4a.  From  it.  we  olitaiu  a  sami)lo  of  OT-ai  points  shown  in  Figure  Mb  (1000  random  points. 
2032  vertex  ])oints.  and  720  boundary  points).  .Mesh  optimization,  with  c,.,,,  =  10“’.  reduces  the 
mesh  down  to  487  v('rti<es  (Figure  1  Ic).  By  setting  r-,.,,,  =  10“'.  we  obtain  a  coar.ser  mesh  of  239 
vertices  ( Figure  1  Id ). 

.■\s  these  examples  illustrate,  basing  mesh  simplification  on  a  measure  of  distance  between  the 
simplified  mesh  and  the  original  has  a  number  of  bemdits; 

•  Vertices  are  dense  in  regions  of  high  (lau.ssian  curvature,  whereas  a  few  large  faces  span  the 
flat  regions. 

•  Long  edges  are  aligned  in  directions  of  low  curvature,  and  the  aspect  ratios  of  the  triangles 
adjust  to  local  curvature. 

•  Edges  and  vertices  of  the  sittiplilit'd  tttesh  are  ])lac('d  in-ar  shar|)  features  of  the  original  mesh. 

5.3  Segmentation 

.Mesh  optimizatiott  (utables  its  to  (huc'cl  sharp  features  iti  the  itttderlyitig  surface.  Fsing  a  simple 
thresholding  method,  the  optimized  tttesh  catt  be  segittettletl  iitto  stttoollt  conti>onett1s.  To  this  end. 
we  build  a  graph  iti  which  the  ttodes  are  the  face.s  of  tttesh.  Two  ttodes  of  this  graph  ;tre  connected  if 
the  two  correspondittg  laces  are  adjacettt  attd  their  dihetiral  attgle  is  stnaller  thait  a  givett  threshold. 
The  connected  cotnjtottettts  of  this  graph  identify  the  desired  stnooth  segntents.  .Vs  ait  example. 
Figure  12e  shows  the  segmeiitatioii  of  t  he  opiimizetl  mesh  into  1  1  components,  .\fter  segmentation, 
vertex  normals  can  be  estimated  from  neighboring  faces  within  t'acli  cotnpotKUtt .  attd  a  smoothly 
shaded  surface  can  be  created  ( Figttre  12f). 

5.4  Parameter  Settings  and  Performance  Statistics 

Table  1  lists  the  specific  parattteter  \alttes  ofc,.,,,  and  h  itsed  to  getterate  the  tneshes  in  the  examples, 
along  with  other  perfortttattce  statistics,  hi  all  I  hes('  examith's.  the  table  tuttry  "vcirud"  refers  to 
a  spring  constant  schedule  of  { L()“-.  10“ ’.  10“ '.  10“'''}.  In  fact,  all  meshes  in  Figure  1  are  also 
created  using  the  same  parameters  (except  that  r,.,,,  was  changt'd  in  two  cases).  I'xecution  times 
were  obtained  on  a  DEC  iiniprocessor  .Alpha  workstation. 
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ral)le  1:  [’('rformaiico  statistics  for  meshes  shown  in  Figure  11. 

6  Related  Work 


Surface  Fitting  There  is  a  large  body  ol' literature  on  lit  ting  embeddings  of  a  rectangular  domain: 
see  Bolle  and  X'emuri  [1]  for  a  review.  Schudy  and  Ballard  [12.  13]  fit  embeddings  of  a  sphere  to 
point  data.  Cio,shtasby  [1]  works  with  embeddings  of  cylinders  and  tori.  Sclaroff  and  Bentland  [14] 
consider  embeddings  of  a  deformed  sui)er(|uadric.  .Miller  el  al.  [.^]  appro.ximate  an  isosurfa.ce  of 
volume  data  by  fitting  a  iik'sIi  liomt'oinorphic  to  a  spher<‘.  While  it  ai)])ears  that  their  method 
could  be  extended  to  finding  isosurla<'es  of  arbil rary  lo|)ological  lyj^e.  it  it  h'ss  obvious  how  it  could 
be  modified  to  handle  point  insttuul  ol  volume  data.  MalK't  [tij  discusses  intei polalion  of  functions 
over  sitnplicial  complexes  of  arbitrary  topological  type. 

Our  method  allows  fitting  of  a  parametric  surface  of  arbit  rary  to|)ological  type  to  a  set  of  three- 
dimensional  points.  In  [2].  we  sketched  an  algorithm  for  fitting  a  mesh  of  fire d  vertex  connectivity 
to  the  data.  The  algorithm  presented  here  is  an  extension  of  this  idea,  in  which  we  also  allow  the 
tiumber  of  vertices  and  their  connectivity  to  vary.  To  the  best  t)f  our  knowledge,  this  has  not  been 
done  before. 


Mesh  Simplification  Two  notal)le  papers  discussing  the  mesh  simplification  probhun  are  Schroeder 
et  al.  [11]  and  Turk  [Ibj. 

The  motivation  of  Schroeder  (U  al.  is  to  simplify  meshes  generated  by  "marching  ctibes"  that 
may  consist  of  more  than  a  million  triangles.  In  their  iterative  aitproach.  the  basic  operation  is 
removal  of  a  vertex  and  re-t  riangulat  ion  of  the  hole  thus  crealc'd.  llu'  criterion  for  ciulex  removal  in 
the  simplest  case  (interior  veutex  not  on  edge  or  corner)  is  the  distance  from  lh('  vertex  to  the  plane 
approximating  its  surrounding  vertices.  It  is  worthwhile  noting  that  this  criterion  only  considers 
deviation  of  the  new  mesh  Irom  the  mesh  created  in  the  |)revious  iteration;  deviatioti  from  the 


original  mesh  does  not.  figure  in  the  strategy. 

Turk’s  goal  is  to  reduce  the  atnount  of  detail  in  a  mesli  wliile  remaining  faithful  to  the  original 
topology  and  geometry.  Ills  Irasic  idea  is  to  distril)ute  points  on  the  existing  mesli  that  are  to 
become  the  new  vertices.  He  then  creates  a  triangulatioii  containing  l)oth  old  and  new  vertices, 
and  finally  removes  tlie  old  vertices.  The  density  of  the  new  vertices  is  chosen  to  be  higher  in  areas 
of  high  curvature. 

The  principal  advantage  of  our  mesh  simplification  method  compared  to  the  tecliniques  men¬ 
tioned  above  is  that  we  cast  mesh  simiilification  as  an  optimization  problem:  we  find  a  new  mesh  of 
lower  complexity  that  is  as  close  as  [rossilrle  to  tlie  original  mesh.  This  is  recognized  as  a  desirable 
property  by  Turk  (.Sectioti  8.  p.  (id):  ".Another  topic  is  finding  measures  of  how  closely  matched  a 
given  re-tiling  is  to  the  original  model.  Can  such  a  (piality  measure  be  used  to  guide  the  re-tiling 
process?".  Optimization  automatically  retains  more  vertices  in  areas  of  high  curvature,  and  leads  to 
faces  that  are  elongated  alotig  directions  of  low  curvature,  another  |)roperty  recognized  as  desirable 
by  Turk. 


7  Summary  and  Future  Work 

We  have  described  an  etiergy  tiiitiiiniztitioti  ap|)ioarh  to  .solving  the  mesh  optimization  problem. 
The  energy  function  we  use  consi.sts  of  t hree  terms:  a  distance  (Uiergy  t hat  meastires  t  he  closeness  of 
fit.  a  representation  energy  that  petializes  meshes  with  a  large  itumber  of  vertices,  and  a  regularizing 
term  that  conceptually  places  spritigs  of  rest  let;gth  zero  on  the  edges  of  the  mesh.  Our  minimization 
algorithm  partitions  the  probleiti  itito  two  nested  subproblems:  an  itiner  cotitinuous  minimizatioti 
and  an  outer  discrete  mitiimiztition.  The  search  s|)ace  consists  of  all  meshes  hoiiK'otnorphic  to  the 
starting  mesh. 

.Mesh  optimization  has  pio\<Mi  elfectiv(>  as  the  s<>cond  jihase  of  our  method  for  surface  recon¬ 
struction  from  unorganized  points,  as  discus.sed  in  [.')].  (  Phase  two  is  r'‘s|)onsible  for  improving  the 
geometric  fit  and  reducing  the  number  of  vert  ices  of  t  he  mesh  produced  in  phase  one.) 

Our  method  lias  also  performed  well  for  mesh  sim|)lilical  ion.  that  is.  the  reduction  of  the  number 
of  vertices  in  a  dens<'  triangubir  mesh.  It  produces  meshes  whose  ('dges  align  themselves  along 
directions  of  low  curvature,  aiul  whose  vertices  concentrate  in  areas  of  high  (laussiau  curvature. 
Because  the  energy  does  not  penalize'  surface's  with  sharp  diheelral  tingles,  the  me'thod  ca'i  recover 
sharp  edges  and  corners. 

•A  number  of  areas  of  future  rese'are  h  still  remain,  including: 

•  Investigate  the  use  of  me)re'  so|)histieate'el  eiptimi/.titiou  nu'thenls.  such  tis  simultiteel  annealing 
for  discrete'  eipt imization  tind  ipitielratie-  me'tlmels  for  uem-linear  letist  seiutires  optimization,  in 
order  to  avoid  unelesirable  loctd  minima  in  the  energy  auel  to  accelerate  convergence. 

•  (lain  more  insight  iiitei  the  ii.se'  eif  the  spring  e-ni'igy  tis  ti  regularizing  term.  es|)e'cially  in  the 
presence  of  a|)preciable  noise. 

•  Improve  thespe'e'el  of  the  algorit  hm  anel  inve'stigate  implementations  on  parallel  architectures. 
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•  Develop  met  hods  lor  liltiiiu,  higher  order  splines  to  more  accurately  and  concisely  model 
curved  surfaces. 

•  Experiment  with  sparse,  non-uniform,  and  noi.sy  data. 

•  Extend  the  current  alf^orithm  to  other  di.stance  measures  such  as  maximum  error  (L^'  norm) 
or  average  error  (i*  norm),  instead  of  the  current  L’  norm. 

A  Mathematical  Appendix 


In  this  appendix,  we  arldress  some  mathematical  issues  alluded  to  in  the  hoily  of  the  paper.  Sec¬ 
tion  .-V.l  contains  a  proof  that  there  is  a  mesh  at  which  the  absolute  minimum  of  the 

energy  function  is  attained.  In  S<'clion  A.'2.  it  is  shown  that  the  sim])licial  operations  introduced 
in  Section  -1.'2  are  comph'te  in  t  lie  scmse  that  they  may  Ix'  used  to  generate  all  triangulations  of 
topological  surfaces.  In  Section  .\.d.  we  prove  tin*  necessity  and  siiinciency  of  the  criteria  for  legality 
of  an  edge  collapse  given  in  Section  1.2 

A.l  Absolute  Minima  for  the  Energy  Functional  Can  be  Attained 

Hecall  that  one  of  the  reasons  for  adding  the  s|)ring  energy  term  to  the  energy  functions 

was  to  insure  the  existence  of  a  mesh  realizing  a  minimum.  In  this  apitendix.  we  prove  that  this  is 
the  case. 

To  .see  what  is  involved,  let  (  >  0  be  the  altsolute  minimum  of  F.  ;md  let  l\\..  ]\.)  be  a  sequence 

of  meshes  with  linn _ x  F{  l\\-.  \\.)  =  < .  Without  the  spring  energy  term,  it  is  itossible  to  construct 

examples  in  which  the  xantex  configuration  \\  Ix'coiik’s  mihoiinded  as  k  ap|)roach('s  inlinity. 

This  cannot  happen  if  the  sprin,s>,  (ux'igy  i<'rm  is  inclmh'd.  We  will  sliow  that  lor  sufficiently 
large  I’  the  energy  fnnclional  /'.(A.  I  )  is  bonmh'd  from  Ixdow  by  a  multiph'  ot  t  lu'  square  ol  the 
Euclidean  norm  of  1’.  ( 'oiise<|U(’iil  !>'.  the  miiiimiz.ing  scMpumci'  (  A'c-  If)  is  cont.'iiiied  in  a  finite  ball 
and  must,  therefore,  contain  a  subsecpience  converging  to  a  mesh  (  A’,„,„.  \ By  construction 
F(K 

T77  ( 71  *  1  nnn  )  —  ^  • 

Theorem  1  Tlicrf  is  a  niish.  { mm)  "'///»  llif  jimixrlii  Ihiil 

FiKA  )  >  /•.(A'., . r„„„) 


for  every  mesh  (  A  .  I  ). 

Remark  1  It  is  imporldnl  to  noh  In  r<  lhat  lln  maij  Or,,,,,,  :  |A',„,,i|  —  R*  nmy  not  lx  an  embeei- 
(linff.  In  (jeneral.  the  minimum  of  F  may  lx  allaiixd  by  a  tlryfnrrale  mai)i)in(j  u'hosr  imaye  may 
not  be  an  embeddal  manifold. 


We  proceed  now  with  the  technicalities  of  I  he  formal  proof. 
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Recall  that  lor  ,■  mesh  (A.l')  and  a  data  set  .V.  //  is  the  miinher  ot  data  points  ,ind  iv  is  the 
number  of  vertices  of  A’.  Assume  that  A’  is  contained  in  a  Itall  of  radius  r  >  0  centered  at  the 

orit^iu.  If  I'  =  (vi.Vi . v,„ )  is  viewed  as  a  vector  in  R’"'  then  its  Euclirlran  norm  on  is 

gi  en  by  the  formula 

l|r|i  = 

where  ||vill  is  the  standard  Euclidean  norm  of  R^. 


Lemma  1  S'k/j/josc  Ihnt  ||\'||  >  4>/m  I'lun 


E{  A',  i  )  >  min 


Proof.  Recall  that 

/■-(  A  .  I  1  =  A.,,„.,,|j(  A .  \  '  4-  £,/;.</(  A  .  ^  )  +  A, ■»/,)(  A  ) . 
Let  b  =  ma.Kj  ||v,||  and  a  =  min,  ||v,|j.  Notice  tlial  b  satisfies  the  ine(|uality 


from  which  the  iiUHiualily  in  the  slateiiient  of  ihe  lemma  I'ollows.  Q.E.D. 

The  next  lemma  is  an  immediate  ronseq nonce  of  I.emma  1  and  the  definition  of  the  energy 
functional. 

Lemma  2  Let  b  =  /■-(  Ao-  Lq)  u'lierr  f  A’u.  V’o)  1“  n  mesh.  Snpitose  that  {K.\']  is  mesh  such 

that  E{h\V)  <  I)  and  let  m  be  the  number  of  vertices  of  A’.  Then 


•  m  <  b/Crep 


Proof,  (of  Theorem  1)  Choose  a  secinence  of  meshes  (A’^  .C/.  )  with  lim  /f(  A \  t )  =  e. 

Let  mu-  he  the  numlter  oi'  \'eiti<'es  of  A’t.  I'iisl  note  that  >in<e  >  c.-fpine.  the  number 

of  vertices  in  h\  Itonnded  l)y  an  integer  .\/.  lint  tluM<'  are  only  a  finite  number  of  simplicial 
complexes  with  at  most  .\l  vertices.  (  onsetinetit Iv.  we  may  choose  a  snlrsetiuence  h with  A’/,.^,  = 
A’/,.^,^i  for  all  j'  and  with  lim,-_.,  lii  ).,)  ~  e.  Let  A'„,,„  denote  such  a  comi>iex. 

By  Lemma  2  the  v(>ctors  all  lie  within  a  irall  of  finite  radius.  By  compactness,  there  is  a 

subsequence  of  \\-  ,  which  converges  to  a  vector  l'In“  mesh  (  A’„„„.  )  satisfies  the  condition 

of  the  theorem.  Q.E.O. 

A. 2  Completeness  of  Operations 

We  want  to  show  that  the  operations  of  edge  split,  edge  swap  and  edge  collapse  form  a  complete 
sot  in  the  sense  that  if  A  and  /,  are  simplicial  snrlac(‘s  with  homeomoi'irhic  i(rpological  realizations 
then  Ij  can  be  obtained  Irom  A  by  a  linit(>  M'<pi('nc('  of  (slue  splits,  edge  swaps  and  edge  collairses. 

.Althotigh  edge  swaps  have  proves!  nsr'fni  iti  onr  optimization  i)rocednre.  the  next  proposition 
shows  that  they  are  not  needed  to  |)rove  completeness. 

Proposition  1  .In  ede/e  siraii  is  i  <iuirale  nt  to  an  eile/e  sjdit  fedloired  bij  an  ide/e  cedlajese . 

Proof.  Let  T\  =  {  c, .  r.}  and  T>  =  {  C| .  cj.  C| }  be  two  I  riangles  of  a  simitlicial  surface  A’  with 

common  edge  e  =  {(q.  cj}.  .\ii  ('dgc'  swap  along  < .  modilie's  A'  to  gi\e  a  new  comitlex  A  '  obtained 
by  replacing  e  with  the  edge  C  =  {c(.C|)  and  T\  and  T<  with  I  ln>  new  triangles  T{  -  {c|.C:j.(q} 
and  =  { i’2 '  i }  riiis  is  (upii valent  to  peuforming  an  edg(>  s|)lit  along  <  followed  by  an  edge 
collapse  (see  Figure  .'').  tj.l',.l). 


Theorem  2  (Completeness)  l,i  I  A  and  L  be  tiro  simidie  ial  surfae'es  siieh  that  |A  |  and  |/.|  are 
hetnie ejinejrphir.  I  hen  I,  is  isonieirphie-  In  ei  sinii>lie  ied  romjde  .v  oblaineel  frenn  A  hi/  ei  Jinile  seeiiience 
of  ede/e  evllapses  and  eele/e  splits. 
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Figure  8:  An  Edge  swaj)  is  o<iuivalcut  to  an  edge  s))lir  I'ollowcd  by  an  edge  collapse. 


To  facilitate  the  jiroof  of  t  lie  t  lieorent.  we  review  of  some  basic  facts  about  simplicial  comple.xes 
(see  [10]).  Recall  that  a  surfarc  K  is  a  liniie  simplicial  compie.x  whose  topological 

realization  |A’|  is  homeomovphic  to  a  contpact  ■i-dimeiisioiial  manil'old.  possibly  with  lioundary.  If 
M  is  a  compact  topological  surface,  a  tridiujuldtion  of  .1/  is  a  simplicial  surface  A',  together  with 
a  homeoniorphism  from  |A’l  to  .1/. 

aabdivisiion  of  A',  written  K'<l\.  is  a  simplicial  comple.x  A  ',  together  with  a  hoineomorphism. 
o  ;  |A’'l  —  lA’].  which  is  alfine  llnoar  on  |n-|  for  each  simph'X  a  6  A'.  Thus,  the  vertices  of  A’'  can 
be  identified  with  [loiiits  of  )A’|  and  A"'  induces  a  triangiilation  of  each  each  2-simple.x  of  |A'|. 

If  A’  and  L  are  two  simjtlicial  surfaces,  a  simpliridi  mt/p.  o  :  A’  —  L  is  a  map  from  the  set  of 

vertice.s  of  A'  to  the  set  of  vertices  of  A',  such  that  if  =  {/■„.  cj . C(,  }  is  a  A-simple.x  of  A'  theit 

o(<t)  =  {o(  C(j).  c5(  ('I ) . o(cr  )}  is  contained  in  a  simple.x  of  /,.  .\  simplicial  map  o  ;  A'  —  A  is 

called  a  i^implicial  isouiariiliisni  if  it  is  a  bij<>ciioii  ;iiid  o~'  is  a  simplicial  map. 

.simplicial  m;ip  o  extends  to  ;i  jiiecewise  linetir  ma|>  |o|  ;  |A  |  —  |A|  by  ilii'  formiil.i 


where  c,  are  the  vmlices-  of  A'  and  6,  are  the  liarycent  tic  coordinates  of  a  point  in  |A'(.  If  o  is  a 
simplicial  isomorphism  then  (o|  is  a  piecewise  linear  liomeomov)>liism. 

We  need  the  following  wi'll  known  iheormn  from  pieci'wisi'  lim'ar  to|)olog\'  [0.  Tln'orem  o.  jtage 
(id) 


Theorem  3  If  1\  find  I.  ar(  hro  sinipliridi  .-iirfdci  irilli  Ihjiik  oniorpliir  lopoUuiicdl  |■((lll:dtions 
then  there  are  fiubdirision.'i  K'  j  A’  and  If  o  /,  and  a  .'•liiipliridl  isoniarphisia  n  ;  A  '  —  1.' . 


Our  i>roof  of  Theorem  2  relies  on  the  ob.servation  that  I  hi‘  operation  of  performing  an  edge 
collapse  can  be  undone  by  a  se(|uence  of  edge  splits  and  edge  collapses. 

^,\s  is  slaiulard.  we  abuse  iiolalioii  sliRlilIv  .iml  iJeiility  verliees,  U-sini|iliees.  aiui  tlieir  images  in  I  he  topological 
realization  of  the  complex.  I  li.il  is.  e  =  (e)  =  |('|. 
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Figure  9:  An  edge  collapse  {(a)  =>  {h))  can  be  reversed  by  a  seciiience  of  edge  splits  {{b)  =>  (r)) 
followed  by  a  secpience  of  edge  collapses  {(c)  =>  (c/)). 

Lemma  3  Siipposf  Ihal  Ihr  siiDpliritil  snrfarr  i  /.s  tihtainnl  from  the  simplirinl  .'surface  A  by  an 
(<lfje  collapse.  Then  K  is  isomorphic  to  a  simplicial  siirfene  obtained  from  L  by  a  si  <ju(  nee  of  edge 
splits  followed  by  a  seepK  m  e  of  edge  collapses. 

Proof.  Suppose  that  L  is  obtained  from  A’  by  collapsiiig  the  edge  e  =  {s.-s'}  which  is  common 
to  the  two  triangles  r  =  {f/.s-.s'}  and  T'  -  {.s..s'./}.  Suppose  t  he  link  of  vertex  s' consists  of  the  set 

of  vertices  {7  =  70- f/i- '/>•  . 'In-I  =  dn+i-  '^}-  I-ot  7',  =  {7,-1  7  =  1 . (»+  H-  Thus.  L 

is  obtained  from  A’  l)y  identifying  s'  with  >  and  |•emoving  the  triangles  T  and  T'  (see  Figure  9  (a) 
and  (b)). 

To  recover  A’  from  /,.  Itegin  by  se<(tientially  performing  edge  splits  along  the  edges  {^.71}. 

{•S.72} . {•s.7,1).  I.abel  the  new  vertices  thus  obtainc'd  >1.  .s_> . s„.  and  let  //  denote  the  final 

simplicial  surface  (see  Figure  9  (c)). 

Lot  A"  denote  llu’  simplicial  surfacr'  obtained  1)>’  se<|nenl  ially  collaitsing  t  lu'  I'dges 

{■S2,.s:5} .  {.s„_i..s„}  and  denoting  l)y  s'  the  single  vertex  to  which  S| . collairse.  It  is  not 

difficult  to  see  that,  with  this  lalieling.  /,"  =  A',  (^ee  Figure  9(d)). 

.\  similar  argniiient  applies  in  the  ca.se  where  e  is  a  boundary  edge  and.  thus,  common  to  only 
one  triangle.  We  leave  a  formal  proof  of  this  case  to  the  ri'adi'r.  (J.F.D. 

Lemma  4  If  h'  is  a  subdirision  of  a  s/mplicial  sarfaee  A  .  then  A  can  be  obtained  from  A'  ria  a 
segucnce  of  edge  collapsis. 

Proof.  Uecausi'  A'  is  a  suixlivision  ol  A  there  is  a  piecc'wise  linear  homeomorphism  o  :  |A-'|  - 
|A'|.  We  will  obtain  A'  from  A"'  in  three  steps: 

Step  1.  For  each  triangle  T  of  A’,  let  \  f-  denote'  tin'  set  of  v(>rtices  of  A'  which  map  under  o 
to  the  interior  of  \r\.  (For  some  triangles.  \  j  may  Ix'  empty.)  Suppose  that  A’l  is  the  simplicial 
surface  obtained  by  identifying  the  xu’itices  of  (>ach  non-i'inply  \  j  with  a  single  \('rf('X  r/-.  This 
can  b('  done  via  a  s('<|uence  of  edge  lollapses  along  <’dges  joining  xcrticc's  forim'd  by  identifying 
vertices  oi  1  7-.  fo  see  that  A’l  is  a  subdivision  of  A’,  we  lu'ed  only  construct  a  jiiecewise  linear 
homeomorphism  o'  :  |A  i|  —  |A  |.  Set  Oi(c)  =  o(e)  if  r  is  not  of  the  form  rj  and  let  o(rj)  be  the 


21 


baryceiiter  of  \T\  for  oacli  r/-.  for  I  7-  ^  0.  It  is  not  <lil!icuU  to  verify  that  Oi  extends  uniquely  to  a 
piecewise  linear  honieoinorithism  Oi  ;  |Ai|  —  |A  |- 

Step  2.  Let  I't,’ denote  the  sot  of  vertices  of  A",  inappin,”;  under  oj  to  t he  interior  oi  an  edge  oi 
|A’|.  Let  I\  2  be  the  siiiiplicial  surface  obtained  Ity  identifying  the  vertices  of  each  non-empty  \  e  to 
a  single  vertex  C£.  Notice  that  A'a  can  be  olrtaincd  from  A'l  by  a  se(]uence  of  edge  collapses  along 
edges  joining  only  pairs  of  vertices  contained  in  a  single  set  of  the  form  1  e-  It  is  not  difficult  to 
show  that  A'2  is  a  subdivision  of  A'  and  that  the  homeomorphism  o-z  :  lANi  —  lA  i  can  be  chosen 
so  that  o"(,ve)  is  the  barycenter  of  the  edge  E  for  each  non-empty  \  e- 

Step  3.  Let  A'.-j  Ire  the  simplicial  complex  obtained  from  AN  by  collapsing  each  edge  joining  a 
vertex  of  the  form  \'e  to  a  vertex  ol  A  .  INicli  vertex  ol  A 7  can  be  identified  with  either  a  \ertex  of 
K  or  with  the  barycenter  of  a  triangle  of  A  and  the  piecewise  linear  map  0\>,  :  A  n  \h  I  induces 
by  this  identification  is  a  homeomor|)hism.  I'lins.  A7  is  a  Milnlivi^ion  ol  A  . 

Step  4.  Finally,  for  each  vertex  of  AN  of  the  form  cj  -  collapse  an  edge  oi  A  t  joining  vj  to  a 

vortex  of  A'.  The  r<'sulting  complex  is  A'.  Q.F.l). 

Proof.(of  Theorem  2)  .  First  note  that  by  Theorem  there  are  subdivisions  A  '  <1  A  and 
V  <  L.  such  that  A''  and  A'  are  isomorphic. 

By  Lemma  -1.  A’  can  be  obtained  from  A'  and  L  Irom  I'  *"•  (inin-  -  ■  edge  collapses. 

But.  by  Lemma  3.  each  edge  collapse  can  be  reversi'd  in  <1  siniuenti  .  .ge  splits  and  edge 

collapses.  Consequently  A  '  can  lie  olitained  from  A'  liy  a  liinte  se<|uence  ol  edge  collapses  and  edge 

splits. 

.Since  there  is  a  finite  se(|uence  of  <'<lge  collap.si's  and  e<lge  splits  translorming  A  into  A  .  which 
is  isomorphic  to  A',  and  there  is  a  finite  sequence  of  edge  splits  and  I'dge  collapses  translorming  L 
into  L.  it  follows  that  there  is  a  finite  se(|ueiice  of  rnlge  splits  and  ('dge  collapses  translorming  A 
into  a  simplicial  complex  which  is  isomorphic  to  A.  (J.l’..l). 

A. 3  Tests  for  Legality  of  an  Edge  Collapse 

The  [mrpose  of  this  appendix  is  to  prove  I  hemaun  I  which  >,ho\\>  itiai  i  he  coiiditicjns  gixen  in 
.Section  4.2.1  are  necessary  and  siidicieni  l<.)r  an  eiige  collapse  10  be  legal. 

Let  A’  l)e  a  simplicial  comi)U>x  whose  topological  reab/.iiioii  iA'j  Is  a  compact  surface  with 
po.s.sibly  non-empty  boundary.  Let  A''  l>e  1  he  simplicial  complex  olttained  by  identifving  t  he  vertices 
i  and  j.  where  {i.j}  is  a  1-simplex  of  A'.  We  say  that  h'  is  obtained  from  A'  by  an  algt  collapse. 
Recall  that  we  retpiire  that  \1\'\  l)e  homeomorphic  to  |A  |.  When  such  is  the  ca.se.  we  say  that  the 
edge  collap.se  {i.j}  —  A  is  Icfjal.  .Vot  all  edge  collapses  are  legal. 

Theorem  4  A<  t  A''  lx  llu  siniplidnl  imopli.v  oblaiii(<l  from  tin  -implicifil  romph  r  A  hp  rollopsiop 
the  edfie  {i.j}.  I  hen  |A’'|  is  honu  omorphir  to  |A|  ij  anil  ooli)  ij  the  folloirinp  cooditions  are  all 
satisfied: 


tl 


/.  A’  has  more  llion  .{  re  liirts  if  oe  ith<  r  {i)  nor  {j}  are  boiinelartj  rertlres.  or  A’  has  more  than 
■]  rertices  if  eillnr  {/}  or  {j}  eire  hounehinj  re  i  tires. 

If  i  anel  j  are  heetli  he)unelanj  rertiees.  {/.j}  is  a  beiiinelarjj  eeleje. 

■?.  For  nil  rertiexs  k  eieljacent  to  bei'h  i  anel  j  (i.e.  {/./.  )  €  A’  emd  {j,h}  €  K ).  {i.j.k}  is  a  face 
of  K. 


We  consider  first  tlie  case  wIkmc  h  has  fewer  than  •')  vertices.  Observe  tlial  no  edge  collapses 
are  legal  if  A’  has  only  d  vertices.  For  each  siinplicial  siirfact'  has  at  least  d  vertices  |  since  it  must 
have  at  least  one  2-siinple.x ).  Since  this  is  true  of  A',  the  complex  h  must  have  at  h'ast  4  vertices. 

Now  suppose  that  A  has  exactly  1  v<'rlices.  I'hen  condition  (d)  is  automatically  satisfied.  If  A 
admits  an  edge  collapst*.  iIkmi  A'  has  d  v<>rtices  and  so  |A’'|  is  homeoinorphic  to  >  ilisk.  But  then 
|A'|  is  also  homeomorphic  to  a  disk  and  so  has  non-cuiipiv  boundary.  .Notice,  that  in  this  case.  A 
may  have  only  two  possible  configurations  feiiher  2  or  .f  fact's),  and  conditions  (  I  )  and  (2)  ensure 
that  only  legal  edge  collapst's  are  allowed. 

We  have  ttow  showti  that  cotidit iotis  (1)  (d)  are  tit'ces.sary  and  siifficietit  in  the  ca.se  where  A’ 
has  fewer  thati  j  vertices,  lleitcefort h.  we  shall  a.ssume  th;it  A'  has  at  least  o  vertices. 

The  proof  of  Theort'iii  I  relies  on  findittg  a  useful  charticterizatioti  of  simjtlicial  surfaces.  To  give 
it.  we  need  a  few  definitiotis.  If  .s  is  a  sitttph'X  of  A'  then  s  denotes  the  interior  of  the  topological 
space  |.s|.  Thus,  if  .s  is  a  2-simplex  theti  s  is  homeomorphic  to  an  open  disk,  if  s  is  a  1-simplex  then 
s  is  homeomorphic  to  ati  opt'ii  ititerval  and  if  .s  i.-.  a  O-simplex  °=  l.sl.  If  L  is  a  subset  of  A’  (not 
tiecessarily  a  subcoitiitlex  )  ihttt  the  /opo/ot/zen/  reeilizeitiem  ol  /,  (alsocalleti  the  mule  rltiinej  toiiolenjical 
space  of  k)  is  the  topological  subspace 


/-I  =  U  ;A  |. 

'el. 


The  stanelarel  eipe  n  elisk  is  the  set  of  points  I)  =  {(.r./y)  |  .i-*  4-  <  1}  C  R"'  and  the  stanelarel 

half-open  elisk  is  the  sul)set  set  /4+  =  {{■>'.  ij)  £  1)  \  ij  >  0}.  The  closures  of  1)  and  in 
are  written  1)  and  l)+.  respectively.  The  stanelarel  ein''^.  written  S.  is  the  boundary  of  D.  The 
half-circle  S+  is  the  intersection  S  n  /2+. 

By  definition,  a  sim|)lici;il  surface  A  is  a  simjzlicial  coniizh'x  with  the  property  thiit.  for  each 
verte.x  r.  there  is  a  homeomorphism  between  (  *  r:  A'|  and  ('ither  D  or  /2+  sending  c  to  the  origin. 
If  V  is  an  interior  vertex  then  |  *  c:  A|  is  hoineomorizhic  to  D  and  if  r  is  a  boundary  vertex  then 
I  ★  V.  A’l  is  homeomorphic  to  /2+.  For  our  purposes,  a  different  (but  ecpovalent )  characterization  is 
needed.  It  is  not  difficult  to  show  that  (  *  c:  A'|  is  homeoinor|)hic  to  D  if  tind  only  if  |link(  c:  A’)|  is 
homeomorphic  to  S  and  that  it  is  hom<'onn)rphic  to  /2+  if  and  oidy  if  |link(c:  A’)|  is  homeomorphic 
to  the  half-circle  .S  +  .  I'liis  leads  to  the  chara(ieri/..itiuu  of  siinplicial  siirfaci's  i;i\(’ii  in  the  next 
lemma. 


Lemma  5  .1  siiiii)li('iiil  (  oiiiijUj-  A’  nf  diiix  n-'^ioii  J  is  a  simj)lit  itil  surface  if  ami  only  if  for  each 
rertex  r  of  K  the  toyoloyical  syaci  \liid\r:  l\  )\  is  home oimjryhir  to  either  the  circle  S  (in  U'hich 
eeise,  r  is  an  interior  rertex)  or  a  half-i  irch  S+  (in  which  case,  r  is  a  boitnelary  rertexj. 


Proof  of  Necessity.  W’e  now  ^liow  iluit  <on(litioiis  (2)  niul  (-'5)  are  necessary  <oiiditions  for 
{i.j}  —  h  lo  be  lef>al. 

Consider  first  condition  (2)  by  sliowing  that  tlie  if  liie  edge  collapse  {i.j}  —  h  is  legal  and  i 
and  j  are  boundary  vertices  then  {Aj}  must  be  a  Iroiindary  edge.  Suppose  that  /  and  j  are  both 
boundary  vertices  and  sujrpose  by  way  of  contradiction  that  {i.j}  is  not  a  Iroundary  edge.  The 
vertices  i  and  j  are  each  incident  iir  two  l)oundary  edges  and  I  lie  boundary  edges  incident  to  i  and 
to  j  are  di.sjoint  (sinci'  thev  necessarily  have  ilisiincl  ('iidpoints).  lint  then  a  total  of  I  boundary 
edges  of  K'  will  be  incident  to  h.  .Sinc<'  h  can  only  he  incitleni  lo  2  boundary  I'dges  \cc  Imve  reached 
a  contradiction. 

To  prove  the  neci'ssily  of  condition  (.{).  .Sup])o.se  that  k  is  a  verie.v  of  A’  such  iliat  {k.i}  and 
{h.j}  are  in  A’  but  {i.j.h}  is  not  a  siiiifile.';  of  A’.  I'liere  are  two  cases  to  consider:  (i)  /,  an  interior 
vertex  and  (ii)  k  a  boundary  verti'X. 

(ij  Suppo.se  that  k  is  an  interior  v(>rtex.  I’hen  |liiik(/.’:  A')|  is  homeomorphic  to  a  circle.  Hence, 
there  are  vertices  a  =  1.2. .. .  //.  a  >  2  such  that 

link(A:  A')  =  {{tn}-{t>:} . {{l>\.  i>>}.  {pi- . ■ 

.Moreover.  +  A  a  simph'x  of  A'  for  a  =  1.2 _ a  (where  +  |  \\'(>  may  arrange 

that  i  =  yi  and  because  {i.j.k}  is  not  a  siiu|)lex  of  A',  j  =  p.,  foi'  sonu'  1  <  n  <  n.  The  existence 
of  such  a  configuration,  howexer.  implii'.-.  that  tin*  space  A'  is  not  a  suilace.  Tor  tin'  link  of  A  in  K’ 
is  the  siinplicial  complex  olilainc'd  from  liiikl/.:  A  )  by  idenillying  the  points  yn  and  p.,.  It  is  easily 
seen  that  the  topological  r<’aii/alion  of  this  i.s  a  "ligitre-';''.  lint  this  contradicts  l.einma  •’)  with  A 
replaced  by  A  '. 

(ii)  Consider  next  the  cast'  where  k  is  a  boundary  xcriex.  Then  then'  are  xi  rtices  a  — 
1.2. ...  II.  II  >  1  such  I  hat 

link(/r:  A')  =  {{p,  }.  {/>.} . {/;„)}  U  {{p\.pj}.  {pi-ps} . {p„^i.p„}}  ■ 

.Moreover.  {/>,,. /a, ^. | . /r }  is  a  simplex  of  A’  for  a  =  1.2 _ a  -  I  and  {/.'./jiI  and  {k.p,,}  are  both 

boundary  edges.  In  this  case.  |link( A:  A  )|  is  homt'otiiorphic  lo  lln*  hall-circle,  .\gain  link(A:  AT 
is  obtained  from  link(A-:  A  )  by  iiieiil ifving  i  ami  j  ami  again  I.emma  .')  applies,  riiiis  jlink( A:  A'')| 
must  be  homeomorphic  to  either  l  lu'  <  ircle  or  the  half-circle. 

If  either  of  i  or  j  is  p,,  for  1  <  u  <  n  then  the  IliitklA-:  A  'll  is  easily  M'en  lo  liotiieomorirhic  to 
neither  the  circle  nor  the  half-circle  (it  is  the  spaced  obtained  by  identifying  an  interior  [loint  of 
t  he  closed  unit  inteiwal  w  it  h  anol  h<’r  |>oitil  <if  i  he  iitimwal  |.  llenci'.  we  may  assume  I  hat  pi  =  /  and 
p„  =  j.  The  union  ol  the  three  edge's  {i.j}.  \j.k}  and  {k.i}  then  forms  a  boundary  comironent  of 
A  .  and  the  elfect  of  the  ('dge  collapse  \i.  j}  —  h  is  lo  rt'tiiovi'  one  of  the  boundary  com|)onents  of 
A.  Hence,  although  it  may  hapiieit  that  A'  is  a  simplicial  surface',  the  two  topological  spaces  |AT 


and  |A’'|  cannot  In*  liomconiorpliic  >in<<*  (il’llK'y  arc  holli  snrraccs)  tlioir  bonndaiics  have  different 
numbers  of  bonndary  components. 

Proof  of  Sufficiency.  .\o\v  suppose  that  the  nnml)er  of  vertices  of  K  is  at  U-ast  5  and  that 
conditions  (2)  and  (d)  are  satislicd.  We  want  to  show  that  the  spaces  |A’|  and  [A'l  are  homeomor- 
phic. 

Observe  that  A'  and  h  '  a‘>;ro('  outside  of  simplicial  neij>.lil)or hoods  of  {/.  _/ }  and  //.  More  precisely, 
if  we  set 

.V  =  */;  A  U  *7:  A'  and  .V'  =  *//.  A  ' 

and 

f  =  An  .V  =  An  A'.  (1) 

I  hen 

A  A’  -  V  =.  A  '  -  V' . 

file  topological  space  |A|  can  thus  be  viewed  in  two  ways: 

•  as  the  space  ol)tained  Iw  |•('moviug  tlie  open  set  |  *  /  U  *j\  from  |A’| 

•  or  as  the  space’  oletained  l).v  removing  the  open  set  |  *  //|  from  |A’'|. 

This  construction  can  lx*  reversed: 

(A'l  =  IA|U|.V|  and  (A-|n(.V|  =  (A'l 

and 

iA"i  =  i/.|u|.V'|  and  |A"ln|.V|  =  Ifj 

(i.e.  |A'|  is  obtained  ley  attacliinn  i.V|  to  |A|  along  the  set  |f '|  and  |A''|  is  obtained  attaching  |.V'| 
to  I  A|  along  \C'\ ). 

Thus.  |A’|  and  |A''|  are  homeomorphic  if  and  only  if  tliere  is  a  liomeomorpliism  between  |A'| 
and  |.V'|  which  is  the’  identity  eui  1( '|.  We  show  ilmi  this  is  in  fact  the  case.  It  is  best  to  consider 
separately  the  three  cases:  (a)  i  and  7  leoth  intc'iior  vertices,  (b)  I’xactly  one  of  i  and  7  are  interior 
vertices  and  (c)  leoth  i  and  7  are  boundary  vertices. 

Case  (a).  Suppose  tliat  /  and  7  are  l)olh  interior  vertices.  I'limi  tliere  are  e.xactly  two  vc’rtices 
l>o  and  r/o  such  that  tlie  siiiiplid’s  \  j.  /./co}  and  are  in  A  .  Because  A’  is  a  simplicial  surface, 

the  links  of  i  and  7  are  circle’s,  ib’iice.  l  h<’r<’  are  verlicc’s 

/A),  in . Pm  =  Vu  Vu-  Vl . <l»  =  /A)  • 

such  that 

{ /.  .  I'  ,+  \}  is  in  A’  for  u  =  0.  1 . ni  -  I 

{ /. f//,. c//,+  i  }  is  in  A’  (dr  />  =  ().! . n  —  1  . 

2.') 


and 


Figure  10:  The  sets  .V  in  the  cases  wliere:  (a)  i  and  j  are  botli  interior  vertices,  (b)  i  is  a  boundary 
verte.x  and  j  is  an  interior  vert(>x.  and  (c)  i  and  j  are  lH)tli  l)oundary  vertices. 

Tlius. 

Iv:  A'l  =  At}l  u  IJ  u 

■  1=0 

and 

_  ,i-\ 

=  M/.J-VolIU  y  A'+i  II  U 

/,=o 

Wo  claim  tliat  the  vertices 

Po-IH . Pm  - 1  •  Vo-  Vi . Va  - 1 

are  di.stiiicl.  To  see  tliis  not<>  first  tliat  a  =  0.  J . in  are  distinct  because  llinkf^’;  A")|  is  a 

circle.  I3y  the  same  reasoning.  <//,.  /;  =  ().! . n  are  distinci.  .\ow  if  //,  =  r/(,  for  some  0  <  n  <  in 

and  0  <  h  <  ii.  then  crindilion  (.()  im|)li<*s  that  {i.j-i>„}  is  in  A’.  Hence,  either  {n.b)  =  (ni.O)  or 
{it.b)  =  (O.n). 

Observe  also  that  in  +  n  >  2.  l'i»r  sn()|)o.se  not  tlnni  {/. /^a.r/o}  and  {J./in-'hi}  are  in  A'  from 
which  it  is  not  diflicnlt  if)  show  ihal  A'  lia.s  only  (  vcnlices  and  is  a  iiU rain'd ron. 

Thus. 

|.'1=  U 

■i=U  '.=() 

is  homeomorirhic  to  D  and  \('\  (the  boundary  of  |.\  |)  is  homeoniorphic  to  .V.  It  is  now  easy  to 
construct  a  homeomorphism  between  |.V|  and  |.V'(  which  is  the  identity  on  |(''|  from  whicli  it  follows 
that  |A’|  and  |A’'|  are  homeoniorphic. 


Case  (b).  .Suppose  that  /  is  a  boundary  vertex  ami  j  is  an  interior  \(’i't(>x.  Then  there  are 
exactly  fv/o  vortices  and  f/n  such  that  {,;. /./A)}  and  \i.j.<in}  are  in  A’.  Because  A’  is  a  simplicial 
surface,  the  link  of  /  is  a  half-circle  and  the  link  of  j  is  a  circh*.  Hence,  there  are  vertices 

'o-  >'] . I'l  =  Pu  •  ami  /;„.  in . p,„  =  f/n  .  and  f/o.  f/i . <i„  . 


such  that 


{ i.  /•.,+  ! }  €  A  for  n  =  I).  I 


t  -  I  . 


and 


{j,Pb^Ph+[}  €  A'  for  />  =  ().! . in  -  1 

{/'.  f/e-  71  + 1 }  €  A'  for  f  =  0.  1 . n  -  1  . 


Consequently. 


and 


m  —  1 


|*j;A'|=  U  l{j.P6.mi}|  U  |{/.j.7o}|  U  |{/.j,;;o}| 

h=0 


1-1 


Tl-1 


1*/:A'|  =  IJ  |{/,r„.r„+i}lU|{Aj,/Jo}lu|{/.j,r/o}|U  U  \{i.flc,<lc+\]\ 


■1=0 


c=0 


We  claim  that  the  vertices 


I'o . I'l-i-Po . Pii-[-‘lo . 7/11  - 1 

are  distinct.  The  reasoniiiii;  is  ilio  saiiK'  as  that  used  in  case  (a):  we  leave  it  to  the  reader  to  fill  in 
the  details.  This  implies  that 


1-1  iii-i  11- 1 

lA-l  =  U  !{'•  I' a  -  /'o+l  }|c  U  U  IJ  l{J.r/c,9c+i}l 

ti=0  h=0  0=0 

is  homeomorphic  to  the  closed  half-disk  D+  and  |C‘j  to  the  half-circle  .S'+.  The  construction  of  a 
homeomorphisin  between  |.V|  and  |A''l  which  is  the  identity  on  |C'|  and  which  sends  i  to  A  is  then 
routine.  Thus,  in  this  case.  too.  tlie  sitaccs  |A'|  and  |A''|  are  liomeomorphic. 

Case  (c)  Reasoning  similar  to  that  of  ca.ses  (a)  and  (1))  sliows  tliat  |.V|  is  liomeomorphic  to  the 
half-disk  Z?+  and  Itf'l  is  homeomor|)liic  to  tlie  half-circle  .V+.  and  that  there  is  a  homeomorphisin 
from  |.V(  to  |.V'|  which  is  the  identity  on  |C|  and  .semis  /  to  h.  Hence.  |A’|  and  |  A''|  are  homeomorphic. 
Details  are  loft  to  the  reader. 


•27 


References 


[1]  Ruud  M.  Bollc  and  Ralia  ('.  Wtiuiri.  Ou  tluoo'diuu'usioual  surface  reconstruction  methods. 

IEEE  PAMI.  .January  1991. 

[2]  T.  DeRose.  II.  Hoppe,  T.  Duchamp.  .1.  .McDonald,  and  \V.  .Stuetzle.  Fitting  of  surfaces  to 
scattered  data.  SPIE.  18;i0;212- >20.  1992. 

[3]  Gene  Golub  and  Charles  Van  Loan.  Matrix  Coinpiitation.s.  .John  Hopkins  l'ni\ersity  Press. 
2nd  edition.  1989. 

[4]  .Vrdeshir  Goshtasby.  .Surface  reconstruction  from  scaitered  measurements.  SPIE.  1830:247- 
2.56.  1992. 

[5]  II.  Hoppe.  T.  DoRo-so.  P.  Duchamp.  .1.  McDotiaUl.  and  W.  Stuetzle.  Surface  reconstruction 
from  unorganized  jtoints.  ( 'oini>iitf  r  draphics  (SKKi'R.X  PH  '02  Prorcfdinf/s).  26(2  );71-78.  .July 
1992. 

[6]  .J.L.  Mallet.  Discrete  smooth  interpolation  in  geometric  modeling.  CAD.  2  K  I ):  178-191.  .April 
1992. 

[7]  Samuel  .Marin  and  J^hilip  Stnit  h.  Parametric  appro.'cimation  of  data  using  ODR  sj)lines.  G.MR 
7057.  General  Motors  Research  l.aborator'tes.  May  1990. 

[8]  .J.V.  Miller.  1).!",.  Hreen.  W.l:].  I.ortui.sen.  R..M.  O'Bara,  and  M..J.  Wozny.  G(‘ometricaJly  de¬ 
formed  models:  .\  metiiod  for  e.\t racting  closed  geometric  models  from  volume  data.  Computer 
Graphics  (SIGGU.APII  '01  /^/WYfd/m/.sy.  25(1  );2 17  226. -Inly  1991. 

[9]  Edwin  E.  .Moi.se.  (iioiintric  Topolatifj  in  Dinx  n.'iions  2  and  -1.  volume  17  of  Graduate  Texts  in 
Mathemtics.  Springer-N’erlag.  .Xew  5'ork.  Hei<lelberg.  Berlin.  1977. 

[10]  C.P.  Rourke  and  B..I.  Sanderson.  I  nt  rodinl  nm  In  pn  <i  irisf  -Inn  ar  Inpe^liM/i/.  \dlume  69  ol  Erejeb- 
nisse  dec  Mathe  niatih  and  direr  (ire  n:</<  hie  h .  Sprim;er-\ ('ilag.  .New  Aork.  Heidelberg.  Berlin. 
1972. 

[11]  William  Schroeder.  .Jotiatlian  Zarge.  and  William  l.orenscm.  Decimation  of  triangle  meshes. 
Computer  Graphics  (SIGGU.APII  '02  PirtceedinejsJ.  2t'){2):{)o-70.  .July  1992. 

[12]  R.  B.  Schudy  ami  I).  H.  Ballard.  Model  di'tection  of  cardiac  chambers  in  ultrasound  images. 
Technical  Report  12.  Computer  Science  Department.  I'niversity  of  Rochester.  1978. 

[13]  R.  B.  Schudy  and  D.  11.  Ballard.  I'owards  an  anatomical  model  of  heait  motion  as  seen  in  4-d 
cardiac  ultrasound  data.  In  Prexe  e  dinejs  nj  the  lilh  (  'nnfe  re  nee  an  Cnmputer  A  pplie'alions  in 
Radiology  and  ( 'ennpule r-.\ieleel  Anedysis  e>J  Reidiedeigiral  Imae/e  s.  1979. 

[14]  Stan  Sclaroffand  .Ale.x  Pentland.  Generalized  iin|)licit  functions  for  computer  graphics.  Com¬ 
puter  Graphirs  (SIGGRAPU  01  Prei<eedingsJ.'2n{A):'2\7  250.  .July  1991. 

[15]  E.  II.  Spaiiier.  Algehraic  Topoleigy.  .McGraw-Hill.  .\’ew  5'ork.  1966. 


[16]  Greg  Turk.  Re-tiling  polygonal  surl'ace.s.  Coinputf  r  Clrapliirs  (SIGGRAPH  '02  Proceedings). 
•26(2):55-64.  July  1992. 

[17]  G.  VVyvill.  (.'.  McPheeters.  anti  13.  W'yvill.  Data  .structures  for  soft  objects.  The  Visual 
Computer.  2(4  ):227-2:34.  August  1986. 


29 


(e)  Optimization  for  fixed  A'n  (k  =  10~^)  (f)  Optimization  with  k  =  10“" 


Figure  11:  Surface  reconstruction  from  simulated  multi-view  range  data. 


(e)  Surface  segmentation  (11  components)  (f)  Smooth  shading  after  segmentation 


Figure  12:  Surface  reconstruction  (continued).  j 
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(e)  Output  of  phase  two  (optimized  mesh)  (f)  Output  of  phase  two  (optimized  mesh) 


Figure  13:  Surface  reconstruction  from  actual  range  data. 


Figure  14:  Mesli  simplification  example 


